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We consider the early time cosmology of f(R) theories in Palatini formalism and study the 
conditions that guarantee the existence of homogeneous and isotropic models that avoid the Big 
Bang singularity. We show that for such models the Big Bang singularity can be replaced by a 
cosmic bounce without violating any energy condition. In fact, the bounce is possible even for 
pressureless dust. We give a characterization of such models and discuss their dynamics in the 
region near the bounce. We also find that power-law lagrangians with a finite number of terms may 
lead to non-singular universes, which contrasts with the infinite-series Palatini f(R) lagrangian that 
one needs to fully capture the effective dynamics of Loop Quantum Cosmology. We argue that these 
models could also avoid the formation of singularities during stellar gravitational collapse. 

PACS numbers: 04.50.Kd, 98.80.-k , 98.80.Qc 



I. INTRODUCTION 



Under the dynamical laws of General Relativity (GR) , 
the standard FLRW cosmology becomes singular at the 
origin of time. The energy density, expansion rate, and 
other geometrical invariants diverge as the volume of 
the Universe goes to zero. The Big Bang singularity 
seems to be an unavoidable and disturbing aspect of 
the currently established cosmological model [lj which, 
to the eyes of many, only a full quantum theory of 
gravity could resolve. Though a completely satisfactory 
quantum theory of gravitation is not yet available, some 
interesting progress has been made in the last years 
in the field of quantum cosmology. The application 
of non-perturbative canonical quantization techniques, 
developed within the framework of Loop Quantum 
Gravity (LQG) j|, to FLRW universes has led to a 
very active research program known as Loop Quantum 
Cosmology (LQC) Q. The chief accomplishment of this 
line of research is the dynamical resolution of the Big 
Bang singularity, which is replaced by a cosmic bounce 
at an energy density of the order of the Planck density. 
Though the dynamics of LQC is governed by difference 
equations, which is reminiscent of the discrete structure 
of spacetime characteristic of LQG, an effective con- 
tinuum description in terms of differential equations is 
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possible. These equations represent modified Friedmann 
dynamics such that the relation between matter/energy 
and curvature is different from that found in GR. As a 
result, at energies close to the Planck scale curvature 
scalars reach an extremum that cannot be dynamically 
surpassed, whereas at low matter /energy densities the 
equations smoothly recover those of GR. Thus, in LQC 
the matter generates curvature but, unlike in GR, it 
cannot break the spacetime. The singularity is, thus, 
resolved by introducing no new degrees of freedom in 
the theory, just by changing the way matter curves the 
geometry. 



These properties of LQC, namely the ability to 
curve spacetime differently than in GR but without 
introducing new degrees of freedom, can also be found 
in modified theories of gravity of the f(R) type in 
Palatini formalism Q (which are radically different from 
the more standard f(R) theories in metric formalism 
[see, for instance, [|[ for a comparison]). This fact was 
first pointed out in Q and was used there to see if a 
Palatini f(R) lagrangian with the same cosmological 
dynamics as LQC could be found [19(. The answer 
was in the affirmative, and a unique solution to the 
problem was found. Though the lagrangian was obtained 
using numerical methods, two approximate analytical 
functions which fit very well the numerical curve at low 
and high curvatures, respectively, were also given. A 
remarkable property of those functions is that they can 
be expressed as infinite power series of the curvature, 
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which puts forward the fact that in non-perturbative 
LQC quantum geometry effects are non-local. 

To the light of the results of [6| one may wonder how 
generic are bouncing cosmologies in Palatini f(R). Can 
we somehow characterize the lagrangians that avoid the 
Big Bang singularity? It is also reasonable to ask if 
such lagrangians are necessarily tied to infinite power 
series, which are associated with the non-local nature 
of curvature found in LQC. These questions and others 
will be explicitly addressed in this paper. We will see 
that the conditions for a lagrangian to generate a cosmic 
bounce are limited and also that finite power series 
lagrangians can produce non-singular cosmologies. In 
fact, we show that the simplest model that one can 
imagine, namely f(R) = R + R 2 / Rp has such a property, 
where Rp is the Planck scale curvature. This model, 
in turn, leads to cosmic bounce even for presureless 
dust, and that makes it particularly interesting for the 
analysis of simple gravitational collapse models. We 
will comment on this aspect and on the extensions with 
respect to the FLRW model needed to correctly address 
the stellar collapse model a la Oppenheimer-Snyder. 



one obtains an algebraic relation between R and T, which 
will be denoted by R — TZ(T) and generalizes the GR 
relation R = —k 2 T. Note that for a given lagrangian one 
may find multiple solutions Hi(T) of {3} whose physical 
validity must be studied model by model. It could also 
happen that (j4|) had no real solutions, which must be 
regarded as unphysical. In physically acceptable models, 
the algebraic relation R = 1Z(T) allows us to see the fp 
term of ([3]) as a function of T. The connection can then 
be solved as the Levi-Civita connection of an auxiliary 
metric hn V conformally related with the physical metric, 
hfiu = fpg^v Using this relation, we can rewrite © as 
follows 



JR 



KfR- f 

2h 



1 

Jr~ 



where the 1Z, fp, and fpp terms on the right hand side 
are functions of T, all covariant derivatives are referred 
to the metric g^ v , and the Einstein tensor on the left 
hand side is also referred to the metric <? M „. 



(5) 



The paper is organized as follows. In section|TT]we dis- 
cuss the FLRW dynamics of Palatini theories and identify 
the conditions for bouncing cosmology to exist. Then we 
consider a simple power-law lagrangian which satisfies 
the bouncing conditions and study its dynamics. In IIVI 
we consider a stellar collapse model following the classical 
treatment of Oppenheimer-Snyder and show that such a 
model is insufficient to guarantee the matching with the 
external Schwarzschild solution at all times. We conclude 
with a summary and some conclusions. 



II. BOUNCING COSMOLOGIES 

The action that defines Palatini f(R) theories is of the 
form 



1 

2^2 



d 4 xy^gf(R) + S m [g^ v ,ip m ] (1) 



Here f(R) is a function of R = g^R^T), with Ra V (T) 
given by R^(T) = -d^ u + d x T^ + T^T^ - lf p r£ A 

where is an independent connection (not the Levi- 
Civita connection of g^ v ). The matter action S m depends 
on the matter fields ip m and the metric g^ v . Variation of 
the action with respect to metric and connection leads to 
the following equations 



t 

where fp denotes df/dR. Using the trace of ([2]), 
Rfp - 2/ = k 2 T , 



(2) 
(3) 

(4) 



When particularized to FLRW cosmologies for matter 
with constant equation of state P = wp, ([5]) becomes 



f R {n 2 p + V/2) 3K 
■'" 



(6) 



(/h + §AiJ 

(2f R + 3A 1 )H - -3H 2 [f R (2 - 3A 2 ) + (2 - 3w)A 1 ] 
+ [tfp(i- w ) + V\-f R (¥p\ (7) 



where H = a/ a, K is the curvature of the spatial homo- 
geneous sections, and we have defined 
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= RfR-f 


(8) 


Ai 


= fpp(a w pRT) 


(9) 


A 2 


fpRR , r> \2 
njpp — J p. 


(10) 


Rt 


K 2 

Rfpp — fp 


(11) 


a w 


= (l + w)(l-3w) 


(12) 



To obtain ^ and ([7]) we have used the conservation law 
p = —3H(1 + w)p and the trace equation ([3]). Note that 
since R = R(p), the functions fp, fpp, V, Ai, A 2 and Rt 
appearing in (|6|) and |7]) are all functions of p. 

One can verify that GR is recovered in the limit of 
linear lagrangian, fpp — > 0. When the non-linearities 
are important, equations and ([7]) represent modified 
Fricdmann equations in which the cosmic expansion is al- 
tered by the new role played by the various p— dependent 
terms. This fact was used in 6] to find an f(R) la- 
grangian that could reproduce the effective dynamics of 
LQC. In the next subsection we summarize those results. 
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Figure 1: Comparison of the ansatz fa given in (|14[> (solid 
line) with the numerical curve obtained in |6| (dashed line). 
The plot represents df /dR versus \R\. 



A. Loop Quantum Cosmology as an f(R) 

The effective dynamics of LQC is characterized by the 
following Friedmann equation (with K = 0) 



3H 2 



K 2 p 



1 



Pcrit 



(13) 



where p cr n ~ 0-^ppianck is the maximum dynamically 
accessible energy density of the theory. An f(R) la- 
grangian reproducing this equation (and also the corre- 
sponding equation for d) can be found by identifying the 
right hand side of (fT3"|) with that of © and replacing p by 
the rule given in ([4|, with T = 2p for a massless scalar. 
One then obtains a second order differential equation for 
the function f(R) which can be numerically solved sub- 
ject to certain boundary conditions, namely agree with 
GR at low curvatures and give the right acceleration a 
at the position of the bounce. The numerical curve was 
very well fitted by the following ansatz (see Fig. 1) 



f(R) 



dR tanh In 

I 103 



R 
l2Rr. 



(14) 



where R c = K 2 p c nt- It is important to note that for 
this isotropic, spatially flat model, the numerical results 
indicate that the bounce occurs when fjj = 0, which here 
corresponds to R = — 12R C . One can readily verify from 
© that fp = is one of the possible conditions for a 
vanishing H when K = 0. Next we discuss the conditions 
for a bounce to exist in the cases with K = —1, 0, +1. 



bounce. We do not consider models with infrared correc- 
tions due to their pathological physical behavior (see the 
discussion in [|| and references therein). For convenience, 
we define A = n 2 p + V/2, which must be positive at low 
curvatures (assuming K/a 2 negligible) because H 2 is al- 
ways positive. In this generic discussion we also assume 
that the term fp + §Ai does not vanish, which would 
imply a divergence of H 2 and, therefore, a Big Bang sin- 
gularity. That possibility will be discussed later, when 
particular models are considered. We also note that the 
particular equation of state w — 1/3 must be excluded 
from the analysis because that particular case implies 
that T = and, therefore, the dynamics of the theory 
is exactly that of GR plus an effective cosmological con- 
stant, which is of no interest in our discussion. 



1. K = case. 

In the spatially flat case there are two possible condi- 
tions for the bounce: 

• Case A: fp = at some R ~ Rp. 

This condition is the one that we find in LQC. From 
H]), we see that if fp — at some high density 
Pb, where the subindex B denotes bounce, then 
f(R B ) = (1 — 3w)k 2 pb/2. This means that at ps 
we find As = 3(1 + w)k 2 ps /4. Since A cannot 
be negative within the interval p £ [0, ps\ (because 
H 2 > 0), then the type A bounce can only occur 
for sources with w > — 1. 

Let us now focus on the region near the bounce. 
Since in the limit fp — > the density reaches its 
maximum value pp, can t> e approximated by 



-(2-3w)H 



2 ' a 2 



(15) 



where a = — B IS a positive constant, 

which follows from the condition for a minimum at 
fp = 0. The solution of (fl~5|) is given by 



H 2 (t) w<i 
H 2 (t) w>i 



2 - 3w 

a 2 
3w-2 



tanh 2 [qV2 - 3w(t - t B )] (16) 



tan 2 [aV3w - 2(t - t B )] (17) 



where tp is an integration constant that sets the 
instant at which the bounce occurs. The corre- 
sponding expansion factor is given by 



B. Characterization of the Bounce. 

We now consider the conditions for getting H = 
which follow from ([5]). We assume that our lagrangian 
is almost linear at low curvatures, i.e., lim/j_»o fp — * 1 
and may depart from unity at high curvatures near the 



a <XL<§ = as(cosh[aV2 - 3w(t - ts)]) 5 ^ (18) 
a(t) w> 2 = a B (cos[oV3u; - 2(t - t B )])^ (19) 

Note that these results are valid for any lagrangian 
f(R) near the bounce caused by the condition fn = 
0. 
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Case B: A = at some R ~ Rp. 

If starting with /r > and A > we assume 
that the bounce occurs when A = 0, which requires 
/(Rb) = RbIrb{^ + 3w)/(3 + 3w), then equation 
dU implies that RbJr b = —3(1 + w)n 2 pB- Re- 
stricting again our attention to the cases w > — 1, 
we see that when A = then RbJr b < 0. If this 
happens keeping fn B > then this condition could 
be physically possible. Otherwise, the condition 
fp = would have been already satisfied and we 
would be in case A. To determine the sign of fp B , 
we focus on the sign of R, which we assume to be 
a monotonous function of p. Assuming that at low 
curvatures the GR relation R = (1 — 3w)k 2 p is a 
good approximation, then i? > if — 1 < to < 1/3 
and R < if w > 1/3 [HJ. Therefore, if w > 1/3 it 
is possible to find lagrangians f(R) for which cos- 
mic bounce occurs before fn = 0. 



2. Cases K 0. 

In this cases, the discussion is not as clean as in the 
K = case, and only a qualitative description is possi- 
ble. For K > 0, we see that the right hand side of ([5]) 
contains a negative term. Therefore, the bounce will oc- 
cur before reaching the conditions fp = or A = of 
the K = case, i.e., at the bounce fn > and A > 0. 
On the contrary, if we choose K < 0, then the curva- 
ture term on the right hand side of ^ is always posi- 
tive and the bounce can only occur if the first term be- 
comes sufficiently negative, which would require either 
f R > 0, A < or fp < 0, A > 0. The latter condition 
is problematic because it implies a change of sign in the 
function fp, which plays a crucial role in the resolution 
of the connection field equation ([3]). In fact, the connec- 
tion is solved as the Levi-Civita connection of an auxi- 
lary metric which is conformally related with the physi- 
cal metric, i.e., h^ v = fpg^y The vanishing of fp then 
raises doubts on the physical validity of such solutions. 
However, we will see later that the numerical evolution 
of such models does not present any pathological behav- 
ior (at least for the models chosen) when the function fp 
crosses the zero. To understand why this solution is also 
physically valid and well behaved, we are forced to inter- 
pret the f(R) theory as a scalar-tensor theory. In fact, it 
is well known that f(R) theories in metric and Palatini 
formalism are dynamically equivalent to the w = and 
w = —3/2 Brans-Dicke theories, respectively, with the 
identifications <j> = fp and V(4>) = Rfp — f . In the Pala- 
tini case, w = —3/2, the scalar field is non-dynamical and 
can be algebraically expressed as a function of the trace 
T, i.e., 4> = <t>(T)- In this formulation of the theory, there 
is no reason, in principle, to forbid a change of sign in 0, 
because one works from the very beginning with g^ and 
there is no need to assume that the conformal transforma- 
tion h^y = (^g^v must always be well defined. One may 
also wonder if the dynamics is well defined at the instant 



when (f> becomes zero. This question is pertinent because 
the field equations can be written as 4>G^ V — K 2 T flL ,+ 
"other terms". It turns out that the "other terms" that 
appear on the right hand side, also contain second deriva- 
tives of the expansion factor which are finite when = 0. 
Therefore, at least for FLRW cosmologies, the dynamics 
is well defined even when = 0. Therefore, these solu- 
tions should also be seen as physically possible. In any 
case, since the current establishment is that K > we 
will not care any more about the possibility K < 0. 



III. TOY MODEL: f(R) = R + 



In this section we consider a simple polynomial model 
with quadratic deviations from GR, f(R) = R + R 2 /Rp. 
Note that Rp is not restricted to be positive. This 
model is particularly interesting because it leads to the 
familiar relation 1Z = —k 2 T, which will make easier the 
discussion and the formulas. We will first focus on the 
case K = to illustrate how the existence of a cosmic 
bounce depends on the choice of equation of state and 
the sign of Rp. Then we consider the case K > 0, which 
can be used as a first approach in the modeling of stellar 
gravitational collapse a la Oppenheimer-Snyder. We will 
show that periodic bounces and expansions occur in this 
model but that a complete description of gravitational 
collapse must go beyond this simple FLRW model for 
the stellar interior (see section HV|) . We will also show 
that cosmic bounce exists for K < 0. 



A. Bouncing Universe with K — 0. 

As we mentioned before, the bounce will occur when- 
ever fp = or A — before the denominator of the right 
hand side of ([6]) can have have the chance to vanish, which 
would imply a singularity. For the model chosen here, (O 
turns into 



H 2 



R 



1 



l-3w TZ 



(1 - 3w) 



1 - (1 + 3w) 



_R_ 

Rp 



(20) 



The condition fp = occurs when TZb = 
translates into n 2 ps = —Rp/2(1 — 3w). 



-Rp/2, which 
From this it 

follows that if Rp > then the bounce can only occur for 
sources with u>>1/3. If w < 1/3 then, we need to flip 
the sign in front of Rp to guarantee the bounce, which is 
equivalent to taking the lagrangian f(R) = R — R 2 /\Rp\. 
In both cases, the denominator does not vanish before 
or at R — TZp if to > — 1. This bouncing mechanism is 
thus possible for sources with w > — 1 (recall that the 
radiation dominated case w — 1/3 is excluded in general) . 



The condition A = 0, occurs when 1Z 
3w), which tranlates into n 2 p^ 



(A) 



-2R P /(l- 
-2R P /(1 - 3w) 2 
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Figure 2: Bouncing cosmologies with K = compared with 
the GR solution (red lines). The continuous black curve rep- 
resents the case w = 0, corresponding to f(R) = R—R 2 /\Rp\. 
The dashed line represents the case to = 1 of the theory 
f(R) = R + R 2 /\Rp\. The initial density is the same for 
the four curves. 



This forces Rp to be negative and implies that only 
f(R) = R — R 2 /\Rp\ could produce a A— induced 
bounce. For this condition to be a valid bouncing 
mechanism, we should have Jr>0 and the denominator 
should not vanish anywhere neither before nor at the 



limiting curvature 1Z 



(A) 



However, it is easy to see 



that fp > can only be satisfied before or at if 
w < — 1 or to > 1/3 while, at the same time, to avoid 

the vanishing of the denominator before or at TZ^ we 
need —1 < w < 1/3. Since these two conditions cannot 
be satisfied together, we conclude that in these models 
there cannot be a A— induced bounce. 



In summary, for the model f(R) = R + R 2 /\Rp\ the 
bounce is possible only for sources satisfying w > 1/3 
whereas for the model f(R) = R — R 2 /\Rp\ the bounce 
is possible only for — 1 < w < 1/3. In both cases, the 
bounce occurs when fp = 0. All other cases lead to 
universes with Big Bang singularities. See Figj2] for 
examples with K = 0. 



B. Cyclic Cosmology with K > 0. 

The condition K > implies that the Universe has a 
finite maximum radius. Since the modified dynamics of 
our model also sets a minimum radius for certain configu- 
rations, the resulting non-singular cosmologies are cyclic. 
We illustrate this cyclic behavior in Fig. [3J where the 
same initial density and models as in the case K = are 
shown. 



Figure 3: Cyclic cosmologies with K > 0. The continu- 
ous black curve represents the case w = 0, corresponding 
to f(R) = R — R 2 /\Rp\. The dashed line represents the case 
it; = 1 of the theory f(R) = R + R 2 /\R P \. The initial density 
is the same as in Fig. [2] 
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Figure 4: Evolution of fn for K = (dashed line), K > 
(dotted line), and K < (solid line). The initial density is 
the same as in Fig. [2] 



C. Bouncing Universe with K < 0. 

The condition K < implies that the bounce occurs 
for negative values of fp. In Fig. 2] we show how the 
sign of fp can change depending on the choice of K. It is 
clear that the numerical curve is smooth and that nothing 
suggests a breakdown of the dynamics when fp = is 
crossed. 



IV. GRAVITATIONAL COLLAPSE 

The last stages of a collapsing massive star were first 
investigated by Oppenheimer and Snyder Q (see also 
fic| for a pedagogical discussion). They simplified the 
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complexity of such astrophysical scenario by assuming 
that the interior of the star could be seen as an homoge- 
neous and isotropic matter distribution. In this manner, 
the interior geometry can be described by means of a 
FLRW metric, whereas the exterior geometry takes the 
form of the Schwarzschild solution. If the collapsing star 
has a finite initial radius, then the interior metric must 
have positive spatial curvature, K > 0. Demanding 
continuity of the metric components across the surface, 
one finds a relation between K , the mass, and the initial 
radius of the star. This simple model showed that for 
an observer comoving with the fluid, the radius of the 
star shrinks to zero and the density goes to infinity in 
a finite proper time. The conclusion was that the fate 
of sufficiently massive stars is the formation of a black 
hole. 

A simplified analysis such as that carried out by 
Oppenheimer and Snyder, in which the matter profile is 
discontinuous at the surface of the star, is not (strictly 
speaking) possible in general for f(R) Palatini theories. 
The reason is that in these theories, the matter must 
satisfy certain differentiability conditions which are not 
needed in the case of GR. This can be understood by 
noting that the field equations ([5]) for the metric in 
Palatini f(R) contain terms of the form V^V^/r, which 
turn into first and second derivatives of the trace T. For 
this reason, a discontinuous matter profile would gener- 
ate divergences [2l| on the right hand side of the metric 
field equations (note that other type of divergences have 
been detected in spherically symmetric stellar solutions 
for some equations of state. Nonetheless, for the models 
considered here they could never be physically realized. 
See the discussion in [ll|, [l2| for details.). Therefore, 
to correctly study a stellar collapse process in Palatini 
f(R) theories, one is (strictly speaking) forced to 
consider an inhomogeneous matter profile that smoothly 
interpolates between the interior distribution of matter 
and the external vacuum. Nonetheless, for the class of 
models that we are considering, these difficulties can be 
relaxed for some time and some qualitative results can 
be extracted. 

Explicit evaluation of the derivatives V^/r, for 
f(R) = R + R 2 1 Rp, shows that those terms are strongly 
suppressed by 1/Rp factors. In fact, in order to get 
iV^/flJ 2 ~ |V M p/p P | 2 , where R P = n 2 p P , on the 
same order of magnitude as n 2 p (the usual Einstein's 
equations dominant term), the matter density profile 
should change by order of the Planck density on scales 
of the order of the Planck length. This means that for 
a realistic initial configuration of macroscopic size, the 
interior can be well approximated as homogeneous and 
isotropic everywhere except at the outermost regions 
of the star, where the density profile must go to zero 
smoothly. If in this inhomogeneous region the length 
scale of variatio n of T is much larger than the Planck 
scale lp ~ l/^/\Rp\, then the radial derivatives V r /i{ 



can be neglected and global internal homogeneity is a 
good approximation. This would justify the use of the 
Oppenheimer-Snyder model to study the first stages 
of the stellar collapse in our theory However, this 
approximate homogeneity cannot last for ever. As the 
star collapses, inhomogeneities must grow faster on 
the outer regions than in the inner ones because radial 
density gradients are larger out there. This should affect 
the rate at which the inner and outer layers of the star 
collapse. One thus expects that the innermost regions, 
where gradients are weaker, should evolve in a way 
similar to that found in the FLRW K > cosmologies 
of FigO This fact suggests that the singularity will be 
avoided. The outermost regions, on the contrary, will 
feel the effect of the density gradients and very little can 
be guessed about their fate. A complete description of 
the matching of these layers of the star with the external 
Schwarzschild solution would thus require detailed 
numerical analysis. The relevant question to answer now 
is whether the surface of the star manages to cross the 
Schwarzschild radius in a finite proper time or not, i.e., 
if an event horizon is formed or not. These issues will be 
studied elsewhere fl3j | . 



V. SUMMARY AND DISCUSSION 

In this work we have shown that Palatini f(R) theo- 
ries have unusual properties that make them specially 
interesting for addressing strong gravity phenomena 
such as the early Universe and stellar collapse processes. 
The simplest lagrangians with high curvature corrections 
that one can consider, namely f(R) = R±R 2 /\Rp\, lead 
to non-singular universes at high curvatures and seem 
compatible with non-singular collapsing stars. Repulsive 
gravitational forces arise in these scenarios induced by 
the time and/or radial derivatives of the matter-energy 
density profile. No new degrees of freedom in the 
gravitational side or in the matter side (exotic sources) 
are needed to get such repulsive gravitational forces. 
This miracle is possible due to the entanglement existing 
between matter and geometry induced by the indepen- 
dent connection used to construct the action. When the 
matter-energy density and pressure reach a certain high- 
energy scale, spacetime is curved differently than in GR 
and singularities are avoided, which contrasts with other 
known mechanisms to avoid cosmic singularities 0. At 
low energies, as expected, the dynamics is identical to 
that of GR. We thus have theories, f(R) = R ± R 2 /R P , 
with the same low energy behavior as GR but without 
its (classical) high-energy drawbacks. 

Quadratic lagrangians of the form R+aR^R^" + PR 2 
have been thoroughly studied in the literature in the 
usual variational metric formalism. It is well known 
from the works of Starobinsky, Schmidt, and Miiller 
[TR [H, [I?], Ell that all (vacuum) solutions of such 



7 



models which recover at late times power-law FLRW 
universes do not avoid the initial Big Bang singularity. 
The only nonsingular solutions of those models describe, 
rather than a bouncing Universe, a universe starting 
from a de Sitter state. When in those models a is set to 
zero, f(R) = R + PR 2 , bouncing solutions can be found, 
though they are incompatible with a power law expand- 
ing Universe at late times (see the discussion in [14| and 
references therein). The Palatini R ± R 2 /\Rp\ models 
discussed here do not suffer from such problems. The 
properties of the generalized models R + aR^R^ + (iR 2 
will be studied elsewhere. 

Our results suggest that Palatini theories might play a 
relevant role in the phenomenology of gravitation at high 
energies. It has already been shown that these theories 
are able to reproduce the dynamics of isotropic LQC. The 
formulation of that theory in terms of holonomics and tri- 
ads and the non-existence of a connection operator seems 
to be related with the fact that an infinite series in R is 
needed to reproduce that dynamics a la Palatini. The 
fact that finite polynomial lagrangians can also lead to 
non-singular universes indicates that non-local dynamics 
is not essential to remove singularities in second-order 
theories. On the other hand, the assumption of indepen- 
dence between metric and connection in the variational 
process is essential to get second order equations for the 



metric. Since full LQG has the same number of dynam- 
ical fields as GR and Palatini f(R), it is thus reasonable 
to expect that effective descriptions of that theory could 
come in the form of Palatini theories. Unfortunately, we 
are still very far from having under control the dynam- 
ical aspects of LQG, which makes impossible the search 
for effective equations. Nonetheless, it would be very 
positive if we could find further relations between Pala- 
tini theories and other aspects of the quantum gravity 
phenomenology, besides that of bouncing cosmologies. A 
better understanding of these theories in scenarios with 
less symmetries and with more general families of Pala- 
tini lagrangians could help us deepen into such aspects. 
That will be the subject of future studies. 
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